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1. Introduction 

Astronomical and cosmological observations have recently shown that Dark Matter 
(DM) is responsible for 23% of the energy budget of the Universe and 83% of its 
mass [l]. The most promising candidate for Dark Matter is the so-called Weakly 
Interacting Massive Particle (WIMP). The existence of WIMPs is independently 
suggested by considerations of Big Bang cosmology and theoretical supersymmetric 
particle phenomenology [21 [3l H] . 

Over the years, many direct detection experiments have been performed to search 
for nuclear recoils due to elastic scattering of WIMPs off the nuclei in the active volume 
of the detector. The main challenge for these experiments is to suppress the backgrounds 
that mimic WIMP-induced nuclear recoils. Today's leading experiments have achieved 
excellent rejection of electromagnetic backgrounds, i.e., photons, electrons and alpha 
particles, that have a distinct signature in the detector. However, there are sources 
of background for which the detector response is nearly identical to that of a WIMP- 
induced recoil, such as the coherent scattering of neutrinos from the sun [5j, or the 
elastic scattering of neutrons produced either by natural radioactivity or by high-energy 
cosmic rays. 

While neutron and neutrino interactions do not limit today's experiments, they are 
expected to become dangerous sources of background when the scale of DM experiments 
grows to fiducial masses of several tons. In traditional counting experiments, the 
presence of such backgrounds could undermine the unambiguous identification of a Dark 
Matter signal because neutrinos are impossible to suppress by shielding and underground 
neutron backgrounds are notoriously difficult to predict [6]. 

An unambiguous positive identification of a Dark Matter signal even in presence 
of unknown amounts of irreducible backgrounds could still be achieved if one could 
correlate the observation of a nuclear recoil in the detector with some unique 
astrophysical signature which no background could mimic. This is the idea that 
motivates directional detection of Dark Matter. 

1.1. The Dark Matter Wind 

The observed rotation curve of our Galaxy suggests that at the galactic radius of the 
sun the galactic potential has a significant contribution from Dark Matter. The Dark 
Matter distribution in our Galaxy, however, is poorly constrained. A commonly used 
DM distribution, the standard dark halo model [7] , assumes a non-rotating, isothermal 
sphere extending out to 50 kpc from the galactic center. The DM velocity is described 
by a Maxwell-Boltzmann distribution with dispersion = 155 km/s. Concentric with 
the DM halo is the galactic disk of luminous ordinary matter, rotating with respect to 
the halo, with an average orbital velocity of about 220 km/s at the radius of the solar 
system. Therefore in this model, an observer on Earth would see a wind of DM particles 
with average velocity of 220 km/s. 

The Dark Matter wind creates two observable effects. The first was pointed out 
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in 1986 by Drukier, Freese, and Spergel [S] who predicted that the Earth's motion 
relative to the galactic halo leads to an annual modulation of the rates of interactions 
observed above a certain threshold in direct detection experiments. In its annual rotation 
around the sun, the Earth's orbital velocity has a component that is anti-parallel to the 
DM wind during the summer, and parallel to it during the winter. As a result, the 
apparent velocity of the DM wind will increase (decrease) by about 10% in summer 
(winter), leading to a corresponding increase (decrease) of the observed rates in DM 
detectors. Unfortunately, this effect is difficult to detect because the seasonal modulation 
is expected to be small (a few %) and very hard to disentangle from other systematic 
effects, such as the seasonal dependence of background rates. These experimental 
difficulties cast a shadow on the recent claimed observation of the yearly asymmetry 
by the DAMA/LIBRA collaboration [9j. 

A larger modulation of the WIMP signal was pointed out by Spergel in 1988. 
The Earth spins around its axis with a period of 24 sidereal hours. Because its rotation 
axis is oriented at 48° with respect to the direction of the DM wind, an observer on 
Earth sees the average direction of the WIMPs change by 96° every 12 sidereal hours. 
This modulation in arrival direction should be resolvable by a Dark Matter directional 
detector, e.g., a detector able to determine the direction of the DM particles. Most 
importantly, no known background is correlated with the direction of the DM wind. 
Therefore, a directional detector could hold the key to the unambiguous observation of 
Dark Matter. 

In addition to background rejection, the determination of the direction of the arrival 
of Dark Matter particles can discriminate [HI [121 [131 [IH [IS] between various DM 
halo distributions including the standard dark halo model, models with streams of 
WIMPs, the Sikivie late-infall halo model [161 Il3 ^S\, and other anisotropic models. 
The discrimination power is further enhanced if a determination of the sense as well 
as the direction of WIMPs is possible [I9]. This capability makes directional detectors 
unique observatories for underground WIMP astronomy. 

1.2. Directional Dark Matter Detection 

When Dark Matter particles interact with regular matter, they scatter elastically off 
the atoms and generate nuclear recoils with typical energies of a few tens of keV, as 
explained in more detail in section O The direction of the recoiling nucleus encodes the 
direction of the incoming DM particle. To observe the daily modulation in the direction 
of the DM wind, an angular resolution of 20-30 degrees in the reconstruction of the 
recoil nucleus is sufficient, because the intrinsic spread in direction of the DM wind is 
~ 45 degrees. Assuming that sub-millimeter tracking resolution can be achieved, the 
length of a recoil track has to be of at least 1-2 mm, which can be obtained by using a 
very dilute g target material. 

An ideal directional detector should provide a 3-D vector reconstruction of the recoil 
track with a spatial resolution of a few hundred microns in each coordinate, and combine 
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a very low energy threshold with an excellent background rejection capability. Such a 
detector would be able to reject isotropy of the recoil direction, and hence identify the 
signature of a WIMP wind, with just a handful of events [I3]. 

More recently. Green and Morgan [19] studied how the number of events necessary 
to detect the WIMP wind depends on the detector performance in terms of energy 
threshold, background rates, 2-D versus 3-D reconstruction of the nuclear recoil, and 
ability to determine the sense of the direction by discriminating between the "head" 
and "tail" of the recoil track. The default configuration used for this study assumes 
a CS2 gaseous TPC running at 0.05 bar using 200 /im pixel readout providing 3-D 
reconstruction of the nuclear recoil and "head-tail" discrimination. The energy threshold 
is assumed to be 20 keV, with perfect background rejection. In such a configuration, 7 
events would be sufficient to establish observation of the WIMP wind at 90% C.L.. In 
presence of background with S/N=l, the number of events necessary to reject isotropy 
would increase by a factor 2. If only 2D reconstruction is available, the required number 
of events doubles compared to the default configuration. "Head-tail" discrimination 
turns out to be the most important capability: if the sense cannot be measured, the 
number of events necessary to observe the effects of the WIMP wind increases by one 
order of magnitude. 

2. Nuclear Recoils in Gaseous Detectors 

To optimize the design of gaseous detectors for directional detection of Dark Matter one 
must be able to calculate the recoil atom energy spectrum expected for a range of WIMP 
parameters and halo models. The detector response in the relevant energy range must 
also be predictable. The response will be governed first and foremost by the track length 
and characteristics (multiple scattering) as a function of recoil atom type and energy. 
Since gas detectors require ionization for detection, design also requires knowledge of 
the ionization yield in gas and its distribution along the track as a function of recoil 
atom type and energy, and possibly electric field. 

The large momentum transfer necessary to produce a detectable recoil in gas implies 
that the scattering atom can be treated as a free particle, making calculations of the 
recoil spectrum essentially independent of whether the target is a solid, liquid, or gas. 

An estimate of the maximum Dark Matter recoil energy for simple halo models is 
given by the kinematically allowed energy transfer from an infinitely heavy halo WIMP 
with velocity equal to the galactic escape speed. This speed is locally about 500-600 
km/sec |20]; WIMPS with higher velocities than this would not be gravitationally bound 
in the halo and would presumably be rare. The corresponding maximum energy transfer 
amounts to < 10 keV/nucleon. The integrated rate will be concentrated at lower energies 
than this, at least in halo models such as the isothermal sphere. For that model, the 
recoil energy (E^) distribution [7] is proportional to exp^—En/Ej), with Ej a constant 
that depends on the target and WIMP masses and the halo model. For a 100 GeV 
WIMP and the isothermal halo model parameters of Ref. [7], Ej/A varies from 1.05 
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to 0.2 keV/nucleon for target mass numbers from 1 to 131. These are very low energy 
particles, well below the Bragg Peak at ~200-800 keV/A. In this regime dE/dx decreases 
with decreasing energy, and the efficiency of ionization is significantly reduced. 

2.1. Lindhard Model for Low-Energy Stopping 

The stopping process for such low energy particles in homoatomicljl substances was 
treated by Lindhard, Scharff, and Schiott [2T], [22] (LSS). This treatment has stood the 
test of time and experiment, making it worthwhile to summarize the results here. 

As is now well-known, the primary energy loss mechanisms for low energy particles 
in matter can be divided into "nuclear stopping", due to atom-atom scattering, and 
"electronic stopping", due to atom-electron scattering. These mechanisms refer only 
to the initial interaction causing the incident particle to lose energy. Nuclear stopping 
eventually contributes to electronic excitations and ionization, and electronic stopping 
eventually contributes to thermal excitations [22J. 

In Ref. [21] the stopping is described using a Thomas-Fermi atom model to obtain 
numerical results for universal stopping-power curves in terms of two variables, the 
scaled energy e = Er/Etf, and the scaled range p = R/Rtf, where Eji and R are 
respectively the energy and the stopping distance of the recoil, and Etf and Rtf are 
scale factor^. 

In Ref. [21] it was shown that nuclear stopping dominates in the energy range where 
most of the rate for Dark Matter detection lies. This can be seen as follows. The scaled 
variables e and p depend algebraically on the atomic numbers and mass numbers of the 
incident and target particles. The scale factor Etf corresponds to 0.45 keV/nucleon for 
homoatomic recoils in Carbon, 1.7 keV/nucleon for Ar in Ar and 6.9 keV/nucleon for 
Xe in Xe. Nuclear stopping ^ was found to be larger than the electronic stopping ^ 
for e < 1.6, which covers the energy range < Ef> < Ej where most of the Dark Matter 
recoil rate can be expected. 

Because of the dominance of nuclear stopping, detectors can be expected to respond 
differently to Dark Matter recoils than to radiations such as x-rays or even a particles, 
for which electronic stopping dominates. Nuclear stopping yields less ionization and 
electronic excitation per unit energy loss than does electronic stopping, implying that 
the W factor, defined as the energy loss required to create one ionization electron, will 
be larger for nuclear recoils. Reference [22] presents calculations of the ultimate energy 
loss partitioning between electronic and atomic motion. Experimenters use empirical 
"quenching factors" to describe the variation of energy per unit of ionization (the "W" 
parameter) compared to that from x-rays. 

I A homoatomic molecular entity is a molecular entity consisting of one or more atoms of the same 
element. 

§ The scale factors are (in cgs-Gaussian units): Etf — ^ZiZT ^^^'^j^^"^ , Rtf — i^^a'^M ^'^m-Mt'' ' Here, 
N= number density of target atoms, subscripts i and T refer to the incident particle and the target 
substance, and a = oq — . '^^^'^ , with oq the Bohr radius. 
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The different microscopic distribution of ionization in tracks dominated by nuclear 
stopping can also lead to unexpected changes in the interactions of ionized and 
electronically excited target atoms (e.g., dimer formation, recombination). Such 
interactions are important for particle identification signatures such as the quantity and 
pulse shape of scintillation light output, the variation of scintillation pulse shape with 
applied electric field, and the field variation of ionization charge collection efficiency. 
Such effects are observed in gases [231 121] , and even more strongly in liquid and solid 
targets [25] . 

Electronic stopping [21] was found to vary as ^ = ky/e with the parameter k 
varying only from 0.13 to 0.17 for homonuclear recoils in A=l to 13 J^. Let us define the 
total stopping as ^ = ^ + ^ and the total scaled range as po = Jq -^Ify. The relatively 
small contribution of electronic stopping and the small variation in k for homoatomic 
recoils, makes the total scaled range for this case depend on the target and projectile 
almost entirely through Etf- 

Predictions for the actual range of homoatomic recoils can be obtained from the 
nearly-universal scaled range curve as follows. Numerically integrating the stopping 
curves of Ref. [21] with k set to 0.15 gives a scaled range curve that fits the surprisingly 
simple expression 

Po = 2.04e + 0.04 (1) 

with accuracy better than 10% for 0.12 < e < 10. According to the formula given 
earlier, the scale factor Rtf lies between 1 and 4 x 10^^ atoms/cm^ for homoatomic 
recoils in targets with 12 < A < 131. Thus the model predicts ranges of several times 
10^^ atoms/cm^ at Eji = Ej. This is of the order of a few mm for a monoatomic gas at 
0.05 bar. As a consequence, tracking devices for Dark Matter detection must provide 
accurate reconstruction of tracks with typical lengths between 1 and a few mm while 
operating at pressures of a small fraction of an atmosphere. 

When comparing LSS predictions with experimental results, two correction factors 
must be considered. First, the widely-used program SRIM [26j produces range-energy 
tables which contain the "projected range", while LSS calculate the path length along 
the track. On the other hand, many older experiments report "extrapolated ranges", 
which are closer in magnitude to the path length than to the "projected range". To 
compare the SRIM tables with LSS, the projected range should be multiplied by a 
factor [21] (1 + where Mt and Mp are the target and projectile masses. This 

correction has generally been applied in the next section, where experimental data are 
discussed. 

In addition, it must be noted that the LSS calculations described above were 
obtained for solids. Therefore, one should consider a gas-solid correction in ranges and 
stopping powers, as discussed by Bohr, Lindhard and Dan ^7\. In condensed phases. 



II rr., 4- ; ■ 0.07932,'° 

i ne parameter fc = — — s4s 

in heavy targets. 



ZiZ2{Ai+A2 
AU2 



1/2 

becomes substantially larger only for light recoils 
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the higher colhsion frequency results in a higher probabihty for stripping of excited 
electrons before they can relax, which leads to a higher energy loss rate than for gases. 
This correction is rather uncertain and has generally not been applied in the following 
section of this paper. 

Finally, numerical calculations to extend the LSS model to the case of targets of 
mixed atomic number are given in Ref. [28] . 

2.2. Experimental Data on Low Energy Stopping in Gases 

The literature of energy loss and stopping of fast particles in matter is vast and 
still growing [29l [30]. However, there is not a lot of experimental data available 
for particle ranges and ionization yields in gas at the very low energies typical of 
Dark Matter recoils, where E/A ~ 1 keV per nucleon. Comprehensive collections of 
citations for all energies are available [26l [31], upon which the widely-used theory- 
guided-fitting computer programs SRIM and MSTAR [31] are based. Several older 
references [321 [331 El] still appear representative of the available direct measurements 
at very low energy. More recent studies [35] provide indirect information based on large 
detector simulations. 

Both references [32] and [33] used accelerated beams of He, N, Ne, Ar and ^^Na, 
^^Ga, and ^^^Au in differentially pumped gas target chambers filled with pure-element 
gases. In [32] the particles were detected with an ionization chamber, while in [33] 
radioactive beams were used. The stopped particles were collected on segmented walls 
of the target chamber and later counted. Typical results were ranges of 2(3.2) x 10^^ 
atoms/cm^ for 26(40) keV Ar+ in Argon. The fit to LSS theory given above predicts 
ranges that are shorter than the experimental results by 10-40%, which is consistent with 
experimental comparisons given by LSS. Accuracy of agreement with the prediction from 
the SRIM code is about the same. As in all other cases discussed below, the direction of 
the deviation from LSS is as expected from the gas-solid effect mentioned in the previous 
section. 

In Ref. [35] nuclear recoils from 252 neutrons were recorded by a Negative Ion 
Time Projection Chamber (NITPC) filled with 40 Torr CS2. The device was simulated 
fitting the observed pulse height and event size distributions. The best fit range curves 
given for C and S recoils in the gas are 10-20% higher at 25-100 keV than LSS predictions 
computed by the present authors by assuming simple additivity of stopping powers for 
the constituent atoms of the polyatomic gas target. 

2.3. Ionization Yields 

Tracking readouts in gas TPC detectors are sensitive only to ionization of the gas. 
As noted above, both nuclear and electronic stopping eventually contribute to both 
electronic excitations (including ionization) and to kinetic energy of target atoms, as 
primary and subsequent generations of collision products interact further with the 
medium. Some guidance useful for design purposes is available from Ref. [22], where 
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the energy cascade was treated numerically using integral equations. In terms of the 
scaled energy e and the electronic stopping coefficient k introduced above, the (scaled) 
energy t] ultimately transferred to electrons was found to be well approximated [36] by 
7] = ^ , with g{e) = e + 3e°'^^ + 0.7e^'^. This function interpolates smoothly from 
?7 = Oate = Oto?7 = e for e — > oo, giving rj = 0.4 at e = 1. In other words, this theory 
predicts only about 40% as much ionization per unit of energy deposited by Dark Matter 
recoils as by low LET radiation such as electrons ejected by x-rays. 

Several direct measurements of total ionization by very low energy particles 
are available in literature. Many of these results are for recoil nuclei from alpha 
decays [STJ [Ml [38]. These ~ 100 keV, A ~ 200 recoils are of interest as backgrounds 
in Dark Matter experiments, but their scaled energy e = 0.07 is below the range of 
interest for most WIMP recoils. Measured ionization yield parameters W were typically 
100-120 eV/ion pair, in good agreement with the approximate formula for rj given 
above. Data more applicable to Dark Matter recoils are given in Refs. |39l [iQl 1311 
Some representative results from these works include [40] W = 91 (65) eV/IP for 25 
(100) keV Ar in Ar, both values about 20% higher than would be predicted by the 
preceding approximate LSS expression. Higher W for gases than for condensed media 
is expected [27] as mentioned above. Ref. [3T] measured total ionization from particles 
with 1 < Z < 22 in methane. While in principle the LSS treatment does not apply 
to heteroatomic gases, using the LSS prescription to predict the W factor for a carbon 
target (rather than methane) yields a value that is 15% lower than the experimental 
results. 

The authors of Ref. [35] also fit their data to derive W-values for C and S recoils. 
Their best-fit values are again 10-25% higher than an LSS-based estimate by the present 
author using additivity. 

To summarize, most of the Dark Matter recoils expected from an isothermal galactic 
halo have very low energies, and therefore nuclear stopping plays an important role. 
The sparse available experimental data on track lengths and ionization yields agrees 
at the ~20% level with simple approximate formulas based on the Lindhard model. 
Without applying any gas-phase correction, LSS-based estimates for range tend to be 
slightly longer than those experimentally measured in gases. The predicted ionization 
parameter W also tends to be slightly lower than the experimental data. This situation 
is adequate for initial design of detectors, but with the present literature base, each 
individual experiment will require its own dedicated calibration measurements. 

3. Considerations for Directional Detector Design 
3.1. Detector Architecture 

From the range-energy discussion in the previous section, we infer that track lengths of 
typical Dark Matter recoils will be only of the order of 0.1 /im in condensed matter, while 
track lengths of up to a few millimeters are expected in gas at a tenth of the atmospheric 
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pressure. Several techniques relevant to direction-sensitive detection using condensed 
matter targets have been reported, including track-etch analysis of ancient mica |43j . 
bolometric detection of surface sputtered atoms , and use of nuclear emulsions |15] . 
The ancient mica etch pit technique was actually used to obtain Dark Matter limits. 
However, recently the focus of directional Dark Matter detection has shifted to low- 
pressure gas targets, and that is the topic of the present review. 

The TPC [I6l HTj is the natural detector architecture for gaseous direction-sensitive 
Dark Matter detectors, and essentially all experiments use this configuration. The active 
target volume contains only the active gas, free of background-producing material. Only 
one wall of the active volume requires a readout system, leading to favorable cost- volume 
scaling. TPCs with nearly 100 m^ of active volume have been built for high energy 
physics, showing the possibility of large active masses. 

3.2. Background Rejection Capabilities 

Gaseous DM detectors have excellent background rejection capability for different 
kinds of backgrounds. First and foremost, direction sensitivity gives gas detectors the 
capability of statistically rejecting neutron and neutrino backgrounds. In addition, 
tracking also leads to extremely effective discrimination against x-ray and 7-ray 
backgrounds [Ml US] . The energy loss rates for recoils discussed in the previous section 
are hundreds of times larger than those of electrons with comparable total energy. The 
resulting much longer electron tracks are easily identified and rejected in any direction- 
sensitive detector. Finally, the measured rejection factors for gamma rays vs. nuclear 
recoils varies between 10^ and 10^ depending on the experiment [50l [35| 151]. 

3. 3. Choice of Pressure 

It can be shown that there is an optimum pressure for operation of any given direction 
sensitive WIMP recoil detector. This optimum pressure depends on the fill gas, the halo 
parameter set and WIMP mass, and the expected track length threshold for direction 
measurement. 

The total sensitive mass, and hence the total number of expected events, 
increases proportionally to the product of the pressure P and the active volume V. 
Equation [1] above shows that the range in atoms/cm^ for WIMP recoils is approximately 
proportional to their energy. Since the corresponding range in cm is inversely 
proportional to the pressure {R oc Er/P), the energy threshold imposed by a particular 
minimum track length Er^min will scale down linearly with decreasing pressure, Er^min 
RminP, where Rmin is the shortest detectable track length. For the exponentially falling 
recoil energy spectrum of the isothermal halo [36] the fraction of recoils above a given 
energy threshold is proportional to exp{—Emin/Eor). Hence the rate of tracks longer 
than the tracking threshold Rmin will scale as iV oc PV exp{—^ RminP) , with ^ a track 
length factor depending on the target gas, WIMP mass, halo model, etc., and the track 
length threshold Rmin depending on the readout technology and the drift distance. This 
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expression has a maximum at Popt = 1/ Rmin] , which shows that the highest event rate 
is obtained by taking advantage of improvement in tracking threshold to run at higher 
target pressure. Operating at this optimum pressure, the track-able event rate still scales 
as PoptV, which increases linearly as the tracking threshold decreases. Achieving the 
shortest possible tracking threshold Rmin is seen to be the key to sensitive experiments 
of this type. 

3.4- Tracking Limit due to Diffusion 

Diffusion of track charge during its drift to the readout plane sets the ultimate limit on 
how short a track can be measured in a TPC Diffusion in gases has a rich phenomenology 
for which only a simplified discussion is given here. More complete discussion with 
references to the literature is given by Rolandi and Blum [52] . 

For low values of electric fields, elementary kinetic theory arguments predict equal 
transverse and longitudinal diffusion to the drift field Ed-, with the rms diffusion spread 
5 given by 



Here k is the Boltzmann constant, T the gas temperature, and L the drift distance. 
No pressure or gas dependence appears in this equation. The diffusion decreases 
inversely as the square root of the applied drift field. Increasing the drift field would 
appear to allow diffusion to be reduced as much as desired, allowing large detectors to 
be built while preserving good tracking resolution. 

However, in reality diffusion is not so easily controlled. The low-field approximation 
given by Equation [2] holds only below a certain maximum drift field value E^"-^, which 
depends on the pressure and target gas. The drift field must not violate the condition 
g^max^ << kT, where the effective mean free path A = 1/fna decreases inversely as the 
pressure. Here a is the average total cross section for scattering of the drifting species 
on the fill gas molecules, n is the number density of molecules, and / is an energy- 
exchange-efficiency factor for the scattering of charge carriers from gas molecules. This 
condition amounts to requiring that the work done by the drift field on a charge carrier 
between collisions and not lost to collisions, must be much smaller than the carrier's 
thermal energy. If this condition is fulfilled it will ensure that the drifting carriers' 
random (thermal) velocity remains consistent with the bulk gas temperature. A larger 
scattering cross section o" or a more effective energy exchange due to strong inelastic 
scattering processes will lead to a shorter effective mean free path and a larger value of 
^max Importantly, iJ™*^^ for electrons in a given gas generally scales inversely as the 
pressure, as would be expected from the presence of the mean free path in the "low 
field" condition. 

If the drift field exceeds E^"-^, the energy gained from the drift field becomes non- 
negligible. The average energy of drifting charge carriers begins to increase appreciably, 




(2) 
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giving them an effective temperature T^fj which can be orders of magnitude larger than 
that of the bulk gas. Under these conditions, the kinetic theory arguments underlying 
equation [2] remain approximately valid if the gas temperature T is replaced by Tefj. 
Diffusion stops dropping with increasing drift field and may rapidly increase in this 
regime, with longitudinal diffusion increasing more rapidly than transverse. 

Values of E^"-^ / P for electrons drifting in various gases and gas mixtures vary from 
~0.1-1 V/cm/Torr at 300 K [531 With drift fields limited to this range and a gas 
pressure of ~ 50 Torr, the rms diffusion for a 1 meter drift distance would be several 
mm, severely degrading the tracking resolution. 

Effects of diffusion can be significantly reduced by drifting negative ions instead of 
electrons [551 EHl ET]. Electronegative vapors have been found which, when mixed into 
detector gases, reversibly capture primary ionization electrons within ~ 100 /im of their 
creation. The resulting negative ions drift to the gain region of the chamber, where 
coUisional processes free the electrons and initiate normal Townsend avalanches |58j . 
Ions have Ejp°^ values corresponding to E/P = 20 V/cm Torr and higher. This is 
because the ions' masses are comparable to the gas molecules, so the energy-exchange- 
efficiency factor / which determines E'^°'^ is much larger than for electrons. Ion-molecule 
scattering cross sections also tend to be larger than electron-molecule cross sections. 
The use of negative ion drift in TPCs would allow sub-millimeter rms diffusion for drift 
distances of 1 meter or larger, although total drift voltage differences in the neighborhood 
of 100 kV would be required. 

The above outline shows that diffusion places serious constraints on the design of 
detectors with large sensitive mass and millimeter track resolution, particularly when 
using a conventional electron drift TPC. 

3.5. Challenges of Directional Detection 

The current limits on spin-independent interactions of WIMPs in the 60 GeV/c^ mass 
range have been set using 300-400 kg-day exposures, for example by the XENONIO [59] 
and CDMS [60] experiments. Next generation non-directional experiments are being 
planned to achieve zero background with hundreds or thousands of times larger 
exposures [61] . 

To be competitive, directional detectors should be able to use comparable 
exposures. However, integrating large exposures is particularly difficult for low-pressure 
gaseous detectors. A fiducial mass of a few tons will be necessary to observe DM-induced 
nuclear recoils for much of the theoretically-favored range of parameter space [1]. This 
mass of low-pressure gas would occupy thousands of cubic meters. It is, therefore, key 
to the success of the directional DM program to develop detectors with a low cost per 
unit volume. Since for standard gaseous detectors the largest expense is represented 
by the cost of the readout electronics, it follows that a low-cost read-out is essential to 
make DM directional detectors financially viable. 
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4. Dark Matter TPC Experiments 

4.I. Early History of Direction-Sensitive WIMP Detectors 

As early as 1990, Gerbier et al. [62] discussed using a hydrogen-filled TPC at 0.02 bar, 
drifting electrons in a 0.1 T magnetic field to detect proton recoils from Dark Matter 
collisions. This proposal was made in the context of the "cosmion", a then-current 
WIMP candidate with very large (10"'^^ cm^) cross section for scattering on protons. 
These authors explicitly considered the directional signature, but they did not publish 
any experimental findings. 

A few years later, the UCSD group led by Masek [63] published results of early trials 
of the first detector system specifically designed for a direction-sensitive Dark Matter 
search. This pioneering work used optical readout of light produced in a parallel plate 
avalanche counter (PPAC) located at the readout plane of a low-pressure TPC. The 
minimum discernible track length was about 5 mm. Electron diffusion at low pressures 
and its importance for the performance of gas detectors was also studied [6l]. This 
early work presaged some of the most recent developments in the field, described in 
section 14.41 

12. DRIFT 

The DRIFT-I collaboration [H] mounted the first underground experiment designed 
for direction sensitive WIMP recoil detection [65]. Re-designed detectors were built 
and further characterization measurements were performed by the DRIFT-II [66] 
collaboration. Both DRIFT detectors were cubical 1 m^ negative-ion-drifting TPCs 
with two back-to-back 0.5 m drift spaces. To minimize material possibly contributing 
radioactive backgrounds, the central drift cathode was designed as a plane of 20 micron 
wires on 2 mm pitch. The endcap MWPCs used 20 /im anode wires on 2 mm-pitch, read 
out with transient digitizers. In DRIFT-II the induced signals on grid wires between 
the MWPC anode and the drift space were also digitized. DRIFT-I had an amplifier- 
and digitizer-per-wire readout, while DRIFT-II signals were cyclically grouped onto a 
small number of amplifiers and digitizers. Both detectors used the negative ion drift 
gas CS2 at nominally 40 Torr, about one eighth of the atmospheric pressure. The 1 m^ 
volume gave approximately 170 grams of target mass per TPC. The CS2 gas fill allowed 
diffusion suppression by running with very high drift fields despite the low pressure. 
DRIFT-II used drift fields up to 624 V/cm (16 V/cm/Torr). 

The detectors were calibrated with alpha particles, ^^Fe x-rays and 252Qf neutrons. 
Alpha particle Bragg peaks and neutron recoil events from sources were quickly seen after 
turn-on of DRIFT-I underground in 2001. Neutron exposures gave energy spectra in 
agreement with simulations when the energy per ion pair W was adjusted in accordance 
with the discussion of ionization yields given above. Simulations of DRIFT-II showed 
that the detector and software analysis chain had about 94% efficiency for detection of 
those 252 neutron recoils producing between 1000 and 6000 primary ion pairs, and 
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a ^°Co gamma-ray rejection ratio better than a few times 10~^ [HZ]. A study of the 
direction sensitivity of DRIFT-II for neutron recoils [68] showed that a statistical signal 
distinguishing the beginning and end of Sulfur recoil tracks ( "head-tail discrimination" ) 
was available, though its energy range and statistical power was limited by the 2 mm 
readout pitch. 

At present two 1 m'^ DRIFT-II modules are operating underground. Backgrounds 
due to radon daughters implanted in the internal surfaces of the detector [67] are under 
study and methods for their mitigation are being developed. The absence of nonzero 
spin nuclides in the CS2 will require a very large increase in target mass or a change of 
gas fill in order to detect WIMPs with this device. 

4.3. Dark Matter Searches Using Micropattern Gas-Gain Devices 

It was shown above that the event rate and therefore the sensitivity of an optimized 
tracking detector improves linearly as the track length threshold gets smaller. In recent 
years there has been widespread development of gas detectors achieving very high spatial 
resolution by using micropatterned gain elements in place of wires. For a recent overview 
of micropattern detector activity, see Ref. [69]. These devices typically have 2-D arrays 
of individual gain elements on a pitch of ~ 0.1 mm. Rows of elements [70] or individual 
gain elements can be read out by suitable arrangements of pickup electrodes separate 
from the gain structures, or by amplifier-per-pixel electronics integrated with the gain 
structure [7T]. Gain-producing structures known as GEM (Gas Electron Multiplier [72] ) 
and MicroMegas (MICRO-MEsh GAseous Structure [73J) have found particularly wide 
application. 

The gas CF4 also figures prominently in recent micropattern Dark Matter search 
proposals. This gas was used for low background work in the MUNU experiment [71] 
and has the advantage of high E^^"-^, allowing relatively low diffusion for electron drift at 
high drift field and reduced pressure [75l [761 ISl] , though it does not approach negative 
ions in this regard. Containing the odd-proton nuclide ^^F is also an advantage since 
it confers sensitivity to purely spin-coupled WIMPs [77j, allowing smaller active mass 
experiments to be competitive. Another attractive feature of CF4 is that its Townsend 
avalanches copiously emit visible and near infrared light [THl [79l [HD] , allowing optical 
readout as in the DMTPC detector discussed in section 14.41 The ultraviolet part of 
the spectrum may also be seen by making use of a wavelength shifter. Finally, CF4 is 
non-flammable and non-toxic, and, therefore, safe to operate underground. 

The NEWAGE project is a current Dark Matter search program led by a 
Kyoto University group. This group has recently published the first limit on Dark 
Matter interactions derived from the absence of a directional modulation during a 
0.15 kg-day exposure ^50]. NEWAGE uses CF4-filled TPCs with a microwell gain 
structure [HH [82l [83]. The detector had an active volume of 23 x 28 x 30 cm^ and 
contained CF4 at 150 Torr. Operation at higher-than-optimal gas pressure was chosen 
to enhance the HV stability of the gain structure. The chamber was read out by a single 
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detector board referred to PIC", preceded by a GEM for extra gas gain. The 

/i-PIC has a micro-well gain structure produced using multi-layer printed circuit board 
technology. It is read out on two orthogonal, 400 micron-pitch arrays of strips. One 
array is connected to the central anode dots of the micro-well gain structure, and the 
other array to the surrounding cathodes. The strip amplifiers and position decoding 
electronics are on-board with the gain structures themselves, using an 8 layer PCB 
structure. 

The detector was calibrated with a 252Qf neutron source. Nuclear recoils were 
detected and compared to a simulation, giving a detection efficiency rising from zero 
at 50 keV to 90% near 250 keV. For comparison, the maximum energy of a ^^F recoil 
from an infinitely heavy WIMP with the galactic escape speed is about 180 keV. The 
measured rejection factor for ^'^''Cs gamma rays was about 10~^. The angular resolution 
was reported as 25° HWHM. Measurement of the forward/backward sense of the tracks 
( "head-tail" discrimination) was not reported. 

Another gaseous Dark Matter search collaboration known as MIMAC [84J is led 
by a group at IPN Grenoble, and has reported work toward an electronically read-out 
direction sensitive detector. They proposed the use of '^He mixtures with isobutane near 
1 bar, and also CF4 gas fills to check the dependence on the atomic number A of any 
candidate Dark Matter signal. The advantages claimed for ^He as a Dark Matter search 
target include nonzero nuclear spin, low mass and hence sensitivity to low WIMP masses, 
and a very low Compton cross section which suppresses backgrounds from gamma rays. 
The characteristic (n,p) capture interaction with slow neutrons gives a strong signature 
for the presence of slow neutrons. The ionization efficiency of ~ 1 keV '^He recoils is 
also expected to be very high, allowing efficient detection of the small energy releases 
expected for this target and for light WIMPs. A micropattern TPC with ~ 350 fim 
anode pitch was proposed to obtain the desired electron rejection factor at a few keV. 
The MIMAC collaboration uses an ion source to generate monoenergetic ^He and F ions 
for measuring the ionization yield in their gas mixtures |85j . 

4.4. DMTPC 

The Dark Matter Time Projection Chamber (DMTPC) collaboration has developed a 
new detector concept [59] that addresses the issue of scalability of directional Dark 
Matter detectors by using optical readout, a potentially very inexpensive readout 
solution. 

The DMTPC detector |86l [87] is a low-pressure TPC filled with CF4 at a nominal 
pressure of 50 torr. The detector is read out by an array of CCD cameras and 
photomultipliers (PMTs) mounted outside the vessel to reduce the amount of radioactive 
material in the active volume. The CCD cameras image the visible and near infrared 
photons that are produced by the avalanche process in the amplification region, 
providing a projection of the 3-D nuclear recoil on the 2-D amplification plane. The 
3-D track length and direction of the recoiling nucleus is reconstructed by combining the 
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measurement of the projection along the amphfication plane (from pattern recognition in 
the CCD) with the projection along the direction of drift, determined from the waveform 
of the signal from the PMTs. The sense of the recoil track is determined by measuring 
dE/dx along the length of the track. The correlation between the energy of the recoil, 
proportional to the number of photons collected in the CCD, and the length of the recoil 
track provides an excellent rejection of all electromagnetic backgrounds. 

Several alternative implementations of the amplification region |5T] were developed. 
In a first design, the amplification was obtained by applying a large potential difference 
(AV = 0.6-1.1 kV) between a copper plate and a conductive woven mesh kept at a 
uniform distance of 0.5 mm. The copper or stainless steel mesh was made of 28 /im 
wire with a pitch of 256 /xm. In a second design the copper plate was replaced with two 
additional woven meshes. This design has the advantage of creating a transparent 
amplification region, which allows a substantial cost reduction since a single CCD 
camera can image tracks originating in two drift regions located on either side of a 
single amplification region. 

The current DMTPC prototype [88] consists of two optically independent regions 
contained in one stainless steel vessel. Each region is a cylinder with 30 cm diameter and 
20 cm height contained inside a field cage. Gas gain is obtained using the mesh-plate 
design described above. The detector is read out by two CCD cameras, each imaging 
one drift region. Two f/1.2 55 mm Nikon photographic lenses focus light onto two 
commercial Apogee U6 CCD cameras equipped with Kodak lOOlE CCD chips. Because 
the total area imaged is 16 x 16 cm^, the detector has an active volume of about 10 
liters. For WIMP-induced nuclear recoils of 50 keV, the energy and angular resolutions 
obtained with the CCD readout were estimated to be ~ 15% and 25°, respectively. This 
apparatus is currently being operated above ground with the goal of characterizing the 
detector response and understanding its backgrounds. A second 10-liter module is being 
constructed for underground operations at the Waste Isolation Pilot Plant (WIPP) in 
New Mexico. 

A 5.5 MeV alpha source from ^^^Am is used to study the gain of the detector 
as a function of the voltage and gas pressure, as well as to measure the resolution as 
a function of the drift distance of the primary electrons to quantify the effect of the 
transverse diffusion. These studies [75l El] show that the transverse diffusion allows 
for a sub-millimeter spatial resolution in the reconstruction of the recoil track for drift 
distances up to 20-25 cm. The gamma ray rejection factor, measured using a ^^^Cs 
source, is better than 2 parts per million [75] . 

The performance of the DMTPC detector in determining the sense and direction of 
nuclear recoils has been evaluated by studying the recoil of fluorine nuclei in interaction 
with low-energy neutrons. The initial measurements were obtained running the chamber 
at 280 Torr and using 14 MeV neutrons from a deuteron-triton generator and a 252Qf 
source. The "head-tail" effect was clearly observed [75l 189] for nuclear recoils with 
energy between 200 and 800 keV. Better sensitivity to lower energy thresholds was 
achieved by using higher gains and lowering the CF4 pressure to 75 torr. These 
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measurements demonstrated [51] "head-tail" discrimination for recoils above 100 keV, 
and reported a good agreement with the predictions of the SRIM [26] simulation. "Head- 
tail" discrimination is expected to extend to recoils above 50 keV when the detector is 
operated at a pressure of 50 torr. To evaluate the event-by-event "head-tail" capability 
of the detector as a function of the energy of the recoil, the DMTPC collaboration 
introduced a quality factor Q{Epi) = e{Eji) x (1 — 2w{Er)Y, where e is the recoil 
reconstruction efficiency and w is the fraction of wrong "head-tail" assignments. The Q 
factor represents the effective fraction of reconstructed recoils with head-tail information, 
and the error on the head-tail asymmetry scales as 1/\J\Q). Early measurements 
demonstrated a Q factor of 20% at 100 keV and 80% at 200 keV [5l]. 

The DMTPC collaboration is currently designing a 1-m^ detector. The apparatus 
consists of a stainless steel vessel of 1.3 m diameter and 1.2 m height. Nine CCD 
cameras and nine PMTs are mounted on each of the top and bottom plates of the 
vessel, separated from the active volume of the detector by an acrylic window. The 
detector consists of two optically separated regions. Each of these regions is equipped 
with a triple-mesh amplification device, located between two symmetric drift regions. 
Each drift region has a diameter of 1.2 m and a height of 25 cm, for a total active volume 
of 1 m^. A field cage made of stainless steel rings keeps the uniformity of the electric 
field within 1% in the fiducial volume. A gas system recirculates and purifies the CF4. 

When operating the detector at a pressure of 50 torr, a 1 m^ module will contain 
250 g of CF4. Assuming a detector threshold of 30 keVee (electron-equivalent energy, 
corresponding to nuclear recoil energy threshold ~ 50 keV), and an overall data-taking 
efficiency of 50%, a one-year underground run will yield an exposure of 45 kg-days. 
Assuming negligible backgrounds, such an exposure will allow the DMTPC collaboration 
to improve the current limits on spin- dependent interactions on protons by about a factor 
of 50 [51]. 

5. Conclusion 

Directional detectors can provide an unambiguous positive observation of Dark Matter 
particles even in presence of insidious backgrounds, such as neutrons or neutrinos. 
Moreover, the dynamics of the galactic Dark Matter halo will be revealed by measuring 
the direction of the incoming WIMPs, opening the path to WIMP astronomy. 

In the past decade, several groups have investigated new ideas to develop directional 
Dark Matter detectors. Low-pressure TPCs are best suited for this purpose if an 
accurate (sub-millimeter) 3-D reconstruction of the nuclear recoil can be achieved. A 
good tracking resolution also allows for an effective rejection of all electromagnetic 
backgrounds, in addition to statistical discrimination against neutrinos and neutrons 
based on the directional signature. The choice of different gaseous targets makes these 
detectors well suited for the study of both spin- dependent (CS2) or spin-independent 
(CF4 and ^He) interactions. 

A vigorous R&D program has explored both electronic and optical readout 
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solutions, demonstrating that both technologies can effectively and efficiently 
reconstruct the energy and vector direction of the nuclear recoils expected from Dark 
Matter interactions. The challenge for the field of directional Dark Matter detection is 
now to develop and deploy very sensitive and yet inexpensive readout solutions, which 
will make large directional detectors financially viable. 
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